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J> ' The interpretation of future precise experiments on atomic parity violation in 

m ' 

terms of parameters of the Standard Model could be hampered by uncertainties in the 
atomic and nuclear structure. While the former can be overcome by measurement in 

m ' 

a series of isotopes, the nuclear structure requires knowledge of the neutron density. 
We use the nuclear Hartree-Fock method, which includes deformation effects, to 
calculate the proton and neutron densities in 125 Cs - 139 Cs. We argue that the good 

a' 

agreement with the experimental charge radii, binding energies, and ground state 
^ ■ spins signifies that the phenomenological nuclear force and the method of calculation 

that we use is adequate. Based on this agreement, and on calculations involving 
different effective interactions, we estimate the uncertainties in the differences of the 
neutron radii 5{r 2 )fq^i and conclude that they cause uncertainties in the ratio of weak 
charges, the quantities determined in the atomic parity nonconservation experiments, 
of less than 10 -3 . Such an uncertainty is smaller than the anticipated experimental 
error. 
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I. INTRODUCTION 



Precision studies of electroweak phenomena provide very important tests of the SU(2)l x 
U(l) Standard Electroweak Model. The measurement of the parity nonconserving (PNC) 
components of the atomic transitions belongs to this class. It offers a unique opportunity for 
testing the electroweak radiative corrections at the one loop level, and, possibly, to search 
for new physics beyond the standard model 

The PNC effects in atoms are caused by the 7, Z° interference in the electron-nucleus 
interaction. The dominant contribution comes from the coupling of the axial electronic 
current to the vector nuclear current. (The interaction of the electronic vector current with 
the nuclear axial current is weaker in heavy atoms, and can be eliminated by summing 
over the PNC effects in the resolved hyperfme components of the atomic transitions. The 
hyperfme dependent effect, which also includes the nuclear anapole moment, is of interest 
in its own right but is not considered hereafter.) Since the vector current is conserved, 
atomic PNC essentially measures the electroweak coupling of the elementary quarks. 

At the present time, PNC measurements in stable 133 Cs atoms have ±2% experimental 
uncertainty ||. (An earlier experiment in Cs was performed by Bouchiat et al. |J; the studies 
of PNC effects in atoms have been reviewed by Commins |7[] and Telegdi ||.) However, 
improvement by an order of magnitude in the experimental accuracy is anticipated and a 
possibility of measuring PNC effects in unstable cesium and francium isotopes has been 
discussed |§ . At this level, two issues must be resolved before an interpretation of the PNC 
data in terms of the fundamental electroweak couplings is possible. The atomic theory, 
even in its presently most sophisticated form |[0|jLl|, introduces about ±1% uncertainty. 



Moreover, the small but non-negligible effects of nuclear size [|T|,[yJ must be addressed. 
This latter problem is the main topic of the present work. 

Atomic PNC is governed by the effective bound electron-nucleus interaction (when taking 
only the part that remains after averaging over the hyperfme components) of the form 
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Hpnc = §fij [~Npn(r) + Z(l - 4 sin 2 9 w )p p {r)} x 4 l5 ^ e d 3 r , (1) 

where the proton and neutron densities p p , n (r) are normalized to unity, and we have assumed 
the Standard Model nucleon couplings 

C lp = 2C lu + C ld = ^(l-Asm 2 6 w ) , (2) 
Cln = Cid + 2C\ U = — - . (3) 



The electron part in Eq. (1) can be parametrized as p!2]JT3 



p 5 (r) ee = C(Z)M(Z, R)f(r) , (4) 

where C(Z) contains all atomic structure effects for a point nucleus, Af is a precisely cal- 
culable normalization factor, and f(r) describes the spatial variation (normalized such that 
/(0) = 1). It is the integrals 



q n , P = J f(r)p PtTl (r)d r (5) 

that determine the effect of the proton and neutron distributions on the PNC observables. 
The for mf actors f(r) can be calculated to the order [Za) 2 for a sharp nuclear surface of 



radius RM.13 



f{r)^l-\{Zaf[{r/Rf -\{r/Rf] . (6) 

For a diffuse nuclear surface numerical evaluation of f(r) is necessary (see below). However, 
the coefficients at (r 2 ) and (r 4 ) remain numerically of the order {Za) 2 and depend only 
weakly on the exact shape of p Pj „(r). In addition, since the electric potential near the 
nucleus is very strong, one can safely neglect atomic binding energies in the evaluation of 
f(r). Below we will separate the effects of the finite nuclear size (i.e., effects related to 
the deviations of q UtP from unity); these terms will be represented by a nuclear structure 
correction to the weak charge. 



Taking the matrix element of Hp^c, one obtains 

(i\H PNC \j) = ^=C(Z)Af[Q w (N, Z) + Q^ C (N, Z)] , (7) 

where Qw{N, Z), the quantity of primary interest from the point of view of testing the 
Standard Model, is the "weak charge." In the Standard Model, with couplings (2) and (3), 
the weak charge is 

Q w = -N + Z(l-Asm 2 6 w ) . (8) 

The nuclear structure correction Q"^ C (N, Z) describes the part of the PNC effect that is 
caused by the finite nuclear size. In the same approximation as Eq. (8) above 

Q™ c = -N(q n - 1) + Z(l - 4 sin 2 6 w )(q p - 1) , (9) 

where q UjP are the integrals of f(r) defined above. (Nuclear structure also affects the normal- 
ization factor TV, which is, however, determined by the known nuclear charge distribution 



mm) 

In a measurement that involves several isotopes of the same element, ratios of the PNC ef- 
fects depend essentially only on the ratio of the weak charges and the corresponding nuclear- 
structure corrections Qw(N, Z) + Q\v c (N, Z). (The dependence Af on the neutron number 
TV will not be considered here.) The ratios of the nuclear-structure corrected weak charges, 
in turn, depend, to a good approximation, only on the differences Aq n of the neutron distri- 
butions in the corresponding isotopes. The uncertainties in these quantities, or equivalently, 
in the differences of the neutron mean square radii 5(A(r 2 ) n,n>), then ultimately limit the 
accuracy with which the fundamental parameters, such as sm 2 9w, can be determined. 

It is the purpose of this work to evaluate quantities q n>p for a number of cesium isotopes, 
which might be used in future high-precision PNC experiments 0. Moreover, we estimate 
the uncertainty in these quantities, respectively in their differences, since they represent the 
ultimate limitations for the interpretation of the PNC measurements. 

In section II we describe the nuclear Hartree-Fock calculations that we performed. In 
section III we compare the calculated binding energies, ground state spins and charge radii 



with the experiment. There we also discuss how corrections for the zero-point vibrational 
motion can be estimated and added. From the spread between the results obtained with two 
different successful effective Skyrme forces, and from the pattern of deviations between the 
calculated and measured isotope shifts in the charge radii, we then estimate the uncertainties 
in the corresponding differences of the neutron radii. Finally, in section IV, we calculate 
the nuclear-structure corrections to the weak charges Q^ C (Z = 55, N = 12 — 84) and their 
uncertainties and discuss the corresponding limiting uncertainties in the determination of 



the fundamental parameters of the Standard Model. (Our notation follows that of Ref. [13 



Others, e.g., Ref. JTOJ] do not explicitly separate the nuclear structure dependent effects. We 
believe that such a separation is very useful, since, as stated above, /(r) in Eq. (||) and 
hence also g np , Eq. (|5|), are essentially independent of atomic structure.) 



II. NUCLEAR HARTREE FOCK CALCULATION 

As demonstrated by numerous calculations, the microscopic description of nuclear ground 
state properties by means of the Hartree-Fock (HF) method with an effective Skyrme force- 
like interaction is remarkably successful The few adjustable parameters in the 
Skyrme force are chosen to fit the various bulk properties (energy per nucleon, compressibil- 
ity modulus, symmetry energy, etc.), and properties of several doubly magic nuclei (binding 
energies, charge radii, etc.) fi~6f . The two most popular sets of Skyrme parameters, namely 
Skyrme III and SkyrmeM* have been successfully employed to describe the properties of 
nuclei in several regions of the periodic table |]T8l , |l9l . Below we show only a few formulae 
essential to the basic understanding of the numerical calculation that we performed; details 
can be found in the quoted references. 

The generalized Skyrme force (including all possible spin-exchange terms and zero-range 
density-dependent interaction) can be written as, 

V 8 = t (l + x P ff )S + + xiP ff )(fc 2 5 + 5k' 2 ) + t 2 (l + x 2 )P CT k ■ 5k' 
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+lhp a 6 + iW{cr x + <x 2 ) • k x 5k' , (10) 
o 

where £0-3,^0-2 and W are the adjustable parameters, and 5 = 5(r — r'). 

Because we are dealing with odd-A nuclei, the unpaired nucleon introduces terms that 
break time-reversal symmetry in the HF functional. When the spin degrees of freedom are 
taken into account, the breaking of time reversal symmetry leads to a rather complicated 
functional f20| , |2T|j . The total energy E, which is minimized in the HF method, can be written 



as a space integral of a local energy density 

E = J H(r)d 3 r , (11) 

with 

h 2 

H(r) = ^r + B lP 2 + B 2 (p 2 n + p\) + B 3 (pr - f) + B 4 (p n r n - f n + p p r p - f p ) 
+B 5P Ap + B 6 (p n A Pn + Pp A Pp ) + B 7 p 2+a + B 8 p a (p 2 n + p\) 
+B 9 (pV ■ J + j -V x p + p n V ■ J n + j n ■ V x Pn + p p V • J p + j p ■ V x p p ) 
+B wP 2 + B n (p 2 n + p 2 p ) + B 12 p a p 2 + B 13 p a (p 2 n + p 2 p ) + E c . (12) 

For complete expressions of the Coulomb energy Ec and the coefficients Bi(i = 1, . . . , 13) 
see Ref. |T[[], where the dependence on Skyrme force parameters in Eq. ([3~Q|) is given. The 
mass densities p T , kinetic density r T , current density j T , spin-orbit density V- J T and vector 
density p T (r = n,p) in Eq. (|12"D can, in turn, be expressed in terms of the single-particle wave 
functions The variation of E with respect to §l{r, a) defines the one-body Hartree-Fock 



hamiltonian h 21 



In the following we will use the mass densities p T , which can be expressed as 

p(r) = E*(^)l 2 • (13) 

k,a 

Here $fc(^*, cr) denotes the component of the kth single-nucleon wave function with spin 
|cr(cr = ±1) along the z direction, and v\ are the BCS occupation factors (see below). The 



expressions for the other densities are again given in Ref. fl21] 



The mean square proton and neutron radii are given by the usual formulae 

J r 2 p T (r)d 3 r . (14) 
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In this work, two discrete symmetries, namely parity and ^-signature, are imposed on 
the wave functions Jl~5] , pT| . The complete description of a wave function requires four real 



functions corresponding to the real, imaginary, spin-up and spin-down parts of HH • 

The numerical approximation to the HF energy E is obtained by a discretization of the 
configuration space on a three-dimensional rectangular mesh. The mesh size Ax is the same 
in the three directions and the abscissae of the mesh points are \{2n + l)Ax. In this work, 
Ax is 0.8 fm, and the mesh size is 16 x 16 x 16. The numerical procedure is described in 



detail in Ref. |L5|. 



Pairing correlations need to be included in a realistic description of medium and heavy 
nuclei. We choose to describe pairing between identical nucleons within the BCS formalism 



using a constant strength seniority force ||15|| . In the usual BCS scheme, the paired states are 
assumed to be the two time-reversed orbitals and Although time reversal symmetry 
is broken in our calculations of odd-A nuclei, the time-reversal breaking terms in the func- 
tional generated by the unpaired odd nucleon are very small compared to the time-reversal 
conserving terms so that the time reversal symmetry is still approximately good. In our cal- 
culation we define the pairing partner of state to be the eigenstate of h whose overlap 
with T$fc is maximal (T is the time-reversal operator). Because the single particle orbital 
occupied by the unpaired nucleon and its signature partner do not contribute to the pairing 
energy, we introduce blocking in our code to prevent these two orbitals from participating 
in pairing and force their BCS occupation numbers to be 1 and 0, respectively. 

As some of the cesium isotopes considered here are deformed, it is very important to 
take the deformation degrees of freedom into account. The method of solving the HF+BCS 
equations by discretization of the wave functions on a rectangular mesh allows any type of 
even multipole deformation. The deformation energy curves are obtained by a constraint 
on the mass quadrupole tensor Qij = (3xiXj — r 2 5ij). The two discrete symmetries of the 



wave functions ensure that the principal axes of inertia lie along the coordinate axes. 
The quadrupole tensor is, therefore, diagonal and its principal values Qi can be expressed 
in terms of two quantities Qo and 7 as 

2 

Qi = <2 cos(7 + i-7r), i = 1,2,3, (15) 
where Q and 7 satisfy the inequalities 

Qo > 0, < 7 < \* ■ (16) 

The values of the three constraints Qi were computed from the desired values of Qo and 7 
and inserted in a quadratic constraint functional added to the variational energy, according 



to the method described in Ref. 221. In the calculations described below, we constrain the 



nuclear shape to be axially symmetric (7 = 0). 



III. COMPARISON WITH EXPERIMENT 

In Fig. 1 we show the potential energy curves for 125 Cs - 139 Cs. According to our calcu- 
lations with Skyrmelll (Skmlll) and SkyrmeM* (SkM*) forces the lighter cesium isotopes 
N < 76 are deformed. For Skmlll such an assignment is able to explain the observed ground 
state spins of | + for N = 70 — 74 and | + for N = 76. For SkM* the mean field proton 
states #7/2 and <i 5 / 2 are interchanged and therefore the ground state spin assignments for the 
deformed cesium isotopes are not correct. (This turns out not to be a very crucial problem.) 
Binding energies and shifts 5r^ n and 5r^ n calculated with the SkM* and Skmlll interactions 
are shown in Tables I and II. The binding energies agree in both cases with the experimental 
values with largest deviation of 4 MeV out of about 1000 MeV of total binding energy. 

The comparison between the measured and calculated isotope shifts is illustrated in 
Figs. 2 and 3 as a series of successively better approximations. First, the crosses, connected 
by dashed lines to guide eyes, show the isotope shifts for spherical nuclei. The agreement 
with experiment is not very good even though the spherical calculation correctly predicts 



S 



that the slope of the dependence Sr 2 (A) is about half of the slope expected from the simple 
relation R = tqA 1 ^ 3 . This means that, on average, the neutron-proton interaction we use 
has the correct magnitude. 

Next, the equilibrium deformation for the lighter cesium isotopes is included (open 
squares), leading to a much better agreement. Further improvement is achieved when the 
effect of zero-point quadrupole vibrational motion is taken into account. It is well known 
that the mean square radius of a vibrating nucleus is increased by 

(r 2 )p=(r 2 )o(l + ^(P 2 )). (17) 

We include this effect of the shape fluctuations using the quantities (f3 2 ) extracted from the 
measured transition matrix elements B(E2, + — > 2 + ) and the relation 

{(3 2 ) = B(E2 1 + -> 2 + )[3ZR 2 /Att]~ 2 . (18) 

We take the average B(E2) of the corresponding Xe and Ba isotopes with neutron numbers 
N = 78 — 84 and correct the radii of 133 Cs- 139 Cs accordingly, as shown in Figs. 2 and 3. 
Thus, further improvement in the comparison with the measured isotope shifts results. (For 
N = 84 the B(E2) values are not known. We use instead the empirical relation between the 
energy of the lowest 2 + state and the deformation parameter B(E2) |[24|| .) This correction 
results in changes in r 2 of 0.2124 fm 2 in 133 Cs, 0.1325 fm 2 in 135 Cs, 0.0724 fm 2 in 137 Cs, and 
0.1263 fm 2 in 139 Cs. 

In a fully consistent calculation, one should make a similar correction for the deformed 
cesium isotopes as well. Since the corresponding B(E2) values for the vibrational states are 
not known, and the corrections are expected to be small, we do not make them. Instead, 
we somewhat arbitrarily assume that the zero-point motion correction is the same as in the 
semimagic 137 Cs. We believe this explains the somewhat poorer agreement in the deformed 
cesium isotopes. 

Even though the quadrupole 2 + states contribute most to the mean square radius via 
Eq. fli~7l), other vibrational states, e.g., the octupole 3~ and the giant resonances, contribute 
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as well; however, all such states not only have smaller collective amplitudes but, even more 
importantly, vary more smoothly with the atomic mass (or neutron number) than the 2 + 
states, and hence their contribution to the shifts 5r 2 should be correspondingly smaller. 

Altogether, the error in the shift 5r 2 is at most 0.2 fm 2 , and appears to be independent of 
the change in the neutron number AN. Thus, for the following considerations we assign an 
uncertainty in the relative value of £r 2 of 0.2 fm 2 . Very little is known experimentally about 
the moments r 4 . Quite conservatively, we assume that the uncertainty in <5r 4 is (r 2 ) x Ar 2 ~ 
5 fm 4 . 

Before turning our attention to the neutron radii, it is worthwhile to make a brief com- 
ment about the comparison with absolute values of (r 2 ) and (r 4 ). Experimentally, muonic 
x-ray energies for the stable 133 Cs have been fitted to the Fermi distribution with the 



halfway radius c = 5.85 fm, surface thickness t = 1.82 fm p5| , p6[ , and (r 2 ) = 23.04 fm . 
Such a Fermi distribution gives (r 4 ) = 673 fm 4 . Our HF calculation corrected for zero- 
point vibrational motion with (/3 2 ) = 0.024, as described above, gives (t^hf — 23.27 fm 2 
for Skmlll and 22.69 fm 2 for SkM* interaction, both quite close to the experimental value. 
The calculated (r 4 ) moments (not corrected for the zero-point motion) are 671(SkmIII) and 
652(SkM*) fm 4 . We see, therefore, that the calculation is quite successful in the absolute 
radii (and even surface thicknesses), in particular for the Skmlll interaction (which gives 
also the correct ground state spin). 

The calculated shifts in the neutron radii 5r 2 and 5r 2 are listed in Tables I (SkM*) and 
II (Skmlll) and the quantities £r 2 corrected for the effect of zero-point vibrational motion 
are displayed in Fig. 4. Several comments about these are in order. First, the slope of the 
dependence of Sr^A) for spherical configurations is correspondingly steeper than the slope 
following from R = r A 1//3 . That is obviously a correct result; the combination of a smaller 
slope in the proton radii and a larger slope in the neutron radii when neutrons are added 
is necessary to maintain on average the R = r^A 1 ^ relation. Second, the HF calculations 
imply that the proton and neutron distributions have essentially identical deformations. This 
agrees with the general conclusion about the isoscalar character of low-frequency collective 
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modes in nuclei (see, e.g., Ref. p7f ). Thus, we accept this result and do not assign any 
additional uncertainty to the possible difference in the deformation of protons and neutrons. 
Finally, for the same reason, we use the same B(E2) values, and the (ft 2 ) extracted from 
them, to correct the neutron radii using Eq. flT7|). Assuming all of the above, we assign 
identical uncertainties to the neutron shifts 5r 2 and the proton shifts 5r 2 , and similarly to 
the fourth moments 5rt „. 

Very little reliable experimental information on neutron distribution in nuclei is available. 
In Ref. iPSfl , data from pionic atoms are analyzed. The corresponding best fit for neutron 
mean square radii agrees very well with the HF results quoted there. The nearest nucleus 
to cesium in Ref. 2£| is 142 Ce. Scaling it with A 2 / 3 one arrives at (r 2 ) = 24.7 fm 2 for 133 Cs, 
somewhat larger than our calculated values 23.7 and 24.0 for SkM* and Skmlll, respectively. 
In Ref. [TIJ, the theoretical neutron density of Brack et al. |29| with (r 2 ) = 23.5 fm 2 , was 
used. That value, presumably obtained by interpolation from the values obtained by the HF 
method using the SkM* interaction, is, not surprisingly, quite close to our calculated values. 
This limited comparison suggests that the absolute radii (r 2 ) have uncertainties of about 
1 fm 2 . The uncertainty in the shifts br 2 n should be substantially smaller, and our estimated 
error of 0.2 fm 2 does not seem unreasonable. 



In Ref. [|13|] the uncertainty in the integrals q ntP was estimated from the spread of the 
calculated values with a wide variety of interactions. Some of the interactions employed in 
13| give better agreement for known quantities (charge radii, binding energies, etc.) than 



others. We chose to use only the two most successful interactions. The spread in the 
calculated shifts Sr 2 n for these two interactions is less than our postulated error of 0.2 fm 2 . 

Pollock et al. Jl3| also argue that the isovector surface term (p p — p n )V 2 (p p — p n ) in 
the Skyrme Lagrangian is poorly determined and may affect the neutron skin significantly, 
without affecting most bulk nuclear properties. We tested this claim by modifying simul- 
taneously the coefficients B 5 — > B 5 (l + x) and B e — ► £> 6 — 2B 5 x in Eq. ([12]). We find that 
when we vary x (i.e., the relative strength of the isovector surface term) from +0.3 to -0.3 
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the proton radius (r 2 ) changes indeed very little (about 0.06 fm 2 ) and the neutron radius 
changes somewhat more (by about 0.1 fm 2 , still less than our estimated error). However, 
the binding energy changes by about 5 MeV, more than the largest discrepancy between 
the theory and experiment. Thus, we do not think that the uncertainty in this particular 
coefficient of the Skyrme force alters our conclusions. 



IV. ESTIMATED UNCERTAINTIES IN PNC EFFECTS 

The nuclear structure effects are governed by the coefficients q UjP , Eq. (|5|), which in turn 
involve integrals of the formfactors /(r), Eq. (|6|). The function f(r) is slowly varying over 
the nuclear volume, and may be accurately approximated by a power series 

/(r) = 1 + f 2 x r 2 + U x r 4 , (19) 

and, therefore, 

q n , P = 1 + / 8 x « p > + f A x « p > . (20) 

For a sharp nuclear surface density distribution the only relevant parameter is the nuclear 
radius R and (r 2n ) = 3/(2n + 3)R 2n . Using the experimental (r 2 ) = 23.04 fm 2 for 133 Cs JZ§, 
we find from Eq. @ 

f(r) = 1 - 2.10 x 10~ 3 r 2 + 1.09 x 10~ 5 r 4 , (21) 

where the distance is measured in fermis. If, instead, we solve numerically the Dirac equation 
for the S1/2 and pi/2 bound electron states in the field of the finite size diffuse surface nucleus, 
we obtain the coefficients /^(/k) of —2.31 x 10~ 3 (1.21 x 10 -5 ) when we use the standard 
surface thickness parameter t = 2.25 fm, and —2.267 x 10 _3 (1.157 x 10~ 5 ) when we use the 
surface thickness adjusted so that the nuclear density parametrized by the two-parameter 
Fermi distribution resembles as closely as possible the Hartree-Fock charge density in 133 Cs. 

The expansion coefficients f2, Ja depend, primarily, on the mean square charge radius. 
To take this dependence into account, we use for 133 Cs the fi and f± above, and for the other 
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isotopes, we use the same surface thickness parameter (t = 1.82) and adjust the halfway 
radius in such a way that the experimental < r 2 > are correctly reproduced. 

It is easy now to evaluate the uncertainty in the factors q n)P given the coefficients f2, f± 
and our estimates of the uncertainties in (r 2 ) and (r 4 ). Substituting the corresponding 
values, we find that the uncertainty is 5q n)P = 4.6xl0 -4 , caused almost entirely by the 
uncertainty in the mean square radii (r 2 ). This uncertainty represents about 1% of the 
deviations of q n , p values from unity. 

Before evaluating the nuclear structure corrections Q^ C (N, Z) we have to consider the 
effect of the intrinsic nucleon structure. Following we use 

« = / d 3 rl(r 2 )Z t , ip , n) f(r)V 2 Pp JQ- n , (22) 

where (r 2 )f nt are the nucleon weak radii, and Q pn are the nucleon weak charges. Neglecting 
the "strangeness radius" of the nucleon, and using the fitted two-parameter Fermi density 
distribution, we find 

q f l = -0.00290, = -0.00102 , (23) 

very close to the sharp nuclear surface values of Pollock et al. fL3"| . The above intrinsic 
nucleon structure corrections are small, but not negligible. More importantly, they are 
independent of the nuclear structure, and cancel out in the differences Aq n>p . 

The quantities 100x(g„-l) and 100x(g p -l) are listed in Table III for all cesium isotopes 
and for the two Skyrme interactions we consider. One can see that they vary by about 
4% for neutrons and are essentially constant for protons when the neutron number increases 
from N = 70 to 84. The variation with N is essentially identical for the two forces, while the 
small difference between the q ntP values calculated with the two forces reflects the difference 
in the absolute values of radii for the two interactions. 

The weak charges Qw(N, Z) and the nuclear structure corrections Qiy c (iV, Z) in Table 
III are radiatively corrected. Thus, instead of the formulae (|j), (||) we use 

Qw(N, Z) = 0.9857 x [-N + Z{1 - 4.012x)] , x = 0.2323 + 0.00365S , (24) 
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following p[ . Here S is the parameter characterizing the isospin-conserving "new" quantum 



loop corrections P(|. Also, 



Q^ C (N,Z) = 0.9857 x [-N{q n - 1) + Z{1 - 4M2x){q p - 1)} . (25) 

These quantities, evaluated for S — 0, are shown in Table III. The assumed uncertainty in 
the shifts of the mean square radii, and consequently in the changes in factors q n)P results 
in the relative uncertainty SQw/Qw of 5xl0 -4 . That uncertainty, therefore, represents the 
"ultimate" nuclear structure limitation on the tests of the Standard Model in the atomic 
PNC experiments involving several isotopes. 

In the atomic PNC experiments involving a single isotope, the uncertainty in the neutron 
mean square radius is larger, and 1 fm 2 appears to be a reasonable choice. Thus, from 
nuclear structure alone, the weak charge in a single isotope has relative uncertainty of 
about 2.5xl0~ 3 , perhaps comparable to the best envisioned measurements, but considerably 
smaller then the present uncertainty associated with the atomic structure. 

Suppose now that in an experiment involving several cesium isotopes one is able to 
determine the ratio 

R(N' N) - Qw(N',Z) + Q^(N',Z) 

i?(iV ' iV) - Q w (N,Z) + Q$r(N,Z) (26) 
with some relative uncertainty 5R/R. To a (reasonable) first approximation 

R(N', N) « Q q^ n [ Z z ] x [1 + qn(N') - q n (N)} . (27) 

Thus, we see that nuclear structure contributes to the uncertainty of R at the level of 
roughly 7xl0~ 4 , where we added the individual errors in quadrature. This uncertainty is 
much smaller than the anticipated experimental error. 

In such a measurement, therefore, the uncertainty in x will be 

5x 5R NN' n 5R 

— w — x w 8 — , (28 

x R ZAN R ' y ' 

(see also [f2"|,|i"5||) where the last factor is evaluated for N',N = 70, 84. The above equation 

illustrates the obvious advantage of using isotopes with large AiV. Also, by performing the 
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measurement with several isotope pairs, one can further decrease the uncertainty 5x. On the 
other hand, the uncertainty in the important parameter S is determined from the relation 
Sx = 0.0036555, and thus 

r „ 5R NN' 

5S m — x — . (29) 

R 0.014ZAiV v ; 

In conclusion, we have evaluated the nuclear structure corrections to the weak charges 
for a series of cesium isotopes, and estimated their uncertainties. We concluded that the 
imperfect knowledge of the neutron distribution in cesium isotopes does not represent in 
the foreseeable future a limitation on the accuracy with which the Standard Model could be 
tested in the atomic PNC experiments. 
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FIGURES 

FIG. 1. The potential energy curves for the isotopes 125 Cs - 139 Cs calculated by the 
Hartree-Fock method using the Skmlll interaction. The notations are shown on the figure. 

FIG. 2. Calculated and experimental isotope shifts 5{r^) in cesium, normalized to the 
semimagic 137 Cs. The SkM* interaction has been used. The correction for zero-point vibrations is 
described in the text. The following notations are used: experimental isotope shift O, spherical HF 
isotope shifts +, HF including equilibrium deformation □, and corrected for zero-point vibrations 
x. 

FIG. 3. Calculated and experimental isotope shifts 5(r 2 ) in cesium. The Skmlll interaction has 
been used. The correction for zero-point vibrations is described in the text. The same notations 
used in Fig. § are used here. 

FIG. 4. Calculated changes in the neutron radii ^(r 2 ) in cesium. The results, corrected for 
zero-point vibrational motion, and calculated with the Skmlll (O) and SkM* (+) interactions, are 
shown. 
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TABLES 

TABLE I. Results of the Hartree-Fock calculations with the SkM* interactions. The exper- 
imental binding energies and isotope shifts 5{r^) are also listed for comparison. (The binding 
energies are in MeV, all radial moments in fm.) The experimental isotope shifts are from Ref. [31], 
normalized to the stable isotope 133 Cs. 
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TABLE II. Results of the Hartree-Fock calculations with the Skmlll interactions. The ex- 
perimental binding energies and isotope shifts 5{r^) are also listed for comparison. (The binding 



energies are in MeV, all radial moments in fm.) The experimental isotope shifts are from Ref. [31] 
normalized to the stable isotope 133 Cs. 
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TABLE III. The radiatively corrected weak charges Qw(N, Z), nuclear structure corrections 
Q^ C (N, Z), and the quantities q n -l, q p -l (the factors (q P:Tl -l) contain the intrinsic nucleon structure 
correction, and are multiplied by 100 for easier display) calculated with the SkM* and Skmlll 
interactions, and with the vibrational corrections described in the text. 
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